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HIGHLIGHTS 


►  The  first  physical  model  for  AEM-DEFCs  is  presented. 

►  The  necessity  of  employing  multi-layer  membrane  model  is  outlined. 

►  Details  of  electroosmotic  and  diffusive  ethanol  crossover  are  discussed. 

►  Fuel  crossover  is  significantly  reduced  in  AEM-DEFCs. 
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A  one-dimensional,  isothermal,  single-phase  model  is  presented  to  investigate  the  mass  transport  in 
a  direct  ethanol  fuel  cell  incorporating  an  alkaline  anion  exchange  membrane.  The  electrochemistry  is 
analytically  solved  and  the  closed-form  solution  is  provided  for  two  limiting  cases  assuming  Tafel 
expressions  for  both  oxygen  reduction  and  ethanol  oxidation.  A  multi-layer  membrane  model  is 
proposed  to  properly  account  for  the  diffusive  and  electroosmotic  transport  of  ethanol  through  the 
membrane.  The  fundamental  differences  in  fuel  crossover  for  positive  and  negative  electroosmotic  drag 
coefficients  are  discussed.  It  is  found  that  ethanol  crossover  is  significantly  reduced  upon  using  an 
alkaline  anion  exchange  membrane  instead  of  a  proton  exchange  membrane,  especially  at  current 
densities  higher  than  500  A  m  2. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

An  extensive,  worldwide  search  for  safe,  renewable  and 
sustainable  sources  of  energy  is  underway.  As  an  alternative  to 
conventional  power  providers  such  as  gas  turbines  and  internal 
combustion  engines  operating  on  a  thermal  cycle,  the  fuel  cell 
technology  has  been  considered  as  a  promising  source  of  clean 
energy.  Fuel  cells  directly  convert  the  chemical  energy  of  fuel  to  the 
electricity  and  are  not  limited  by  the  maximum  Carnot  efficiency. 
Hydrogen  is  the  most  commonly  used  fuel  in  fuel  cells,  especially  in 
a  proton  exchange  membrane  type.  Issues  regarding  storing, 
handling,  transporting  and  production  of  hydrogen  led  us  to 
consider  energy-dense  and  hydrogen-rich  liquid  alcohols  such  as 
methanol  and  ethanol. 

Methanol  is  the  most  attractive  alternative  due  to  its  large 
hydrogen  to  carbon  ratio  in  the  molecular  structure.  During  the 
past  decade,  extensive  research  efforts  have  been  focused  on  the 
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study  of  the  direct  methanol  fuel  cell  (DMFC)  [1,2],  In  addition  to 
the  obstacle  of  methanol  crossover  [3,4]  and  the  sluggish  kinetics  of 
methanol  oxidation  [5],  the  toxicity  of  methanol  presents  a  major 
disadvantage.  Ethanol,  on  the  other  hand,  offers  a  higher  specific 
energy  (30  MJ  kg-1)  than  methanol  (19.7  MJ  kg'1)  and  is  safe  and 
renewable.  Ethanol  is  the  mass  produced  by  fermentation  of  sugar- 
containing  biomass  such  as  sugar  cane,  corn  and  wheat. 

Early  efforts  to  use  ethanol  were  based  on  DMFC  technology  and 
using  Nation  as  a  membrane  in  an  acidic  medium.  It  demonstrated 
a  poor  performance  compared  to  methanol,  nearly  1/7  of  maximum 
power  density  of  methanol  at  room  temperature  [6],  Recently,  the 
novel  idea  of  changing  the  reacting  medium  from  an  acid-base  to 
alkali-base  has  made  a  breakthrough  improvement  in  employing 
ethanol  as  a  sustainable  energy  source  [6],  The  alkaline  medium  is 
provided  by  employing  an  anion  exchange  membrane  (AEM). 
Utilizing  an  alkaline  medium  versus  an  acidic  medium  offers 
several  benefits  including:  (i)  improved  oxygen  reduction  and 
ethanol  oxidation  which  allows  the  use  of  non-noble  catalysts,  (ii) 
better  corrosion  resistance,  and  (iff)  negative  electro-osmotic  drag, 
mitigated  fuel  crossover  and  resulting  less  potential  water  flooding 
in  the  cathode. 
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Nomenclature 

Greek 

Transfer  coefficient 

AEM 

Anion  Exchange  Membrane 

7 

Reaction  order 

AEM-DEFC  Direct  ethanol  fuel  cell  using  an  AEM 

e 

Porosity  of  porous  medium 

C 

Concentration,  molm-3 

V 

Overpotential,  V 

CL 

Catalyst  layer 

P 

Density,  kg  m  3 

DEFC 

Direct  ethanol  Fuel  Cell 

a 

Conductivity,  0  1  m  1 

D 

DL 

Diffusivity,  m2s-1 

Diffusion  layer 

<l> 

Phase  potential,  V 

DMFC 

Direct  Methanol  Fuel  Cell 

Superscripts 

Erev 

Reversible  cell  voltage,  V 

CL 

Catalyst 

F 

Faraday  constant,  96,485.3,  A  s  mol-1 

eff 

effective 

G 

H 

Gibbs  energy,  J 

Enthalpy,  J 

ref 

Reference  value 

I 

Current  density,  A  m  2 

Subscripts 

b 

Parasitic  current,  A  m-2 

a 

Anode 

jref 

Jofl2 

Reduction  exchange  current  density,  A  m  3 

AEM 

amb 

Anion  Exchange  Membrane 

Ambient 

Jot 

Oxidation  exchange  current  density,  A  m  3 

c 

Cathode 

L 

Length,  m 

car 

Carbon 

n 

Number  of  computational  node  (layer) 

CL 

Catalyst 

n" 

Molar  flux,  mol  m~2s-1 

cr 

Crossover 

nd 

Electroosmotic  coefficient 

d 

Drag 

R 

Reaction  rate,  A  m-3 

D 

Diffusive 

Ru 

Universal  gas  constant,  8.3145  J  mol_1K-1 

e 

Ethanol  or  electron 

Rcon 

Contact  resistance,  Q  m2 

E 

Electroosmosis 

Rr 

Internal  resistance,  0  m2 

el-ch 

Electrochemical 

s 

Entropy,  J  K-1 

h2o 

Water 

s 

Liquid  saturation 

KOH 

Potassium  hydroxyl 

T 

Temperature,  K 

lim 

Limiting  value 

u 

Velocity,  m  s-1 

mem 

Membrane  phase 

u 

Thermodynamic  equilibrium  potential,  V 

o 

Standard  condition  of  T  =  298  K  and  P  —  1  atm 

V 

Voltage,  V 

ox 

Oxidation 

X 

Coordinate  for  mass  transport  in  the  cell,  m,  or  mole 

02 

Oxygen 

fraction 

OH 

Hydroxyl  ion 

X' 

Coordinate  for  charge  transport  in  the  anode  catalyst 

P 

Parasitic 

layer,  m 

red 

Reduction 

X 

Coordinate  for  charge  transport  in  the  cathode  catalyst 

res 

Reservoir 

layer,  m 

0 

x'  =  0  or  x"  =  0 

During  the  past  few  years,  there  have  been  increasing  research 
interests  in  various  aspects  of  a  direct  ethanol  fuel  cell  using  an 
AEM,  known  as  an  AEM-DEFC.  Table  1  presents  a  summary  of  the 
experimental  studies  conducted  on  the  AEM-DEFC  performance 
and  design. 

The  table  is  broken  down  into  active  and  passive  AEM-DEFCs.  In 
an  active  system,  fuel  and  oxidant  are  forced  through  the  channels, 
whereas  in  a  passive  system  there  are  no  moving  parts.  There  are,  in 
general,  two  types  of  AEMs:  (i)  AEMs  with  a  quaternary  ammonium 
functional  group  and  (ii)  polybenzimidazole  (PB1)  doped  in  strong 
bases.  The  former  is  usually  known  as  the  name  of  its  major 
supplier,  Tokuyama,  and  the  latter  benefits  from  the  higher  ion 
conductivity  (up  to  9  Q-1m-1 )  [19],  Tokuyama  membrane  (denoted 
by  T  in  the  second  column  of  Table  1 ),  Pd-based  catalyst  for  the 
anode  and  Fe— Co  based  catalysts  for  the  cathode  are  most 
commonly  used.  Adding  a  hydroxyl  source  to  the  inlet  fuel  is 
necessary  to  increase  the  ion  conductivity  of  the  AEM  and  to 
improve  the  kinetics  of  ethanol  oxidation  reaction  (EOR).  Both  KOH 
and  NaOH  may  be  added  to  the  anode  side.  In  all  the  studies  con¬ 
ducted,  it  was  preferable  to  use  KOH  rather  than  NaOH  as  the 
hydroxyl  source  in  the  anode  channel,  because  KOH  solution  has 
a  higher  ionic  conductivity  than  a  NaOH  solution  under  the  same 
conditions.  More  importantly,  both  K2CO3  and  CH3COOK  have 


much  higher  solubility  in  the  water  than  Na2C03  and  G^COONa 
[20],  The  highest  reported  power  density  for  an  AEM-DEFC  is  170 
mW-cm-2  at  80  °C  using  pure  oxygen  in  the  cathode  for  the  active 
mode  [11]  and  65  mW-cm-2  at  80  °C  using  air  in  the  cathode  for 
the  passive  mode  [18], 

To  the  authors’  knowledge,  there  is  no  reported  physical  model  for 
transport  phenomena  in  an  AEM-DEFC.  From  a  practical  perspective, 
thorough  knowledge  of  mass  transport,  which  is  difficult  to  obtain 
through  the  experimental  measurements,  is  of  enormous  importance 
for  the  further  development  of  this  type  of  fuel  cell.  The  objective  of 
this  study  is  to  provide  detailed  insight  concerning  the  mass  transport 
in  an  AEM-DEFC  incorporating  the  negative  electroosmotic  drag.  A 
multi-layer  membrane  model  is  proposed  to  investigate  the  funda¬ 
mental  influences  of  using  an  alkaline  membrane  in  reducing  ethanol 
crossover.  The  electrochemistry  is  analytically  solved,  and  the  closed- 
form  solution  is  provided  for  two  limiting  cases.  Ethanol  crossover  is 
substantially  reduced  when  an  AEM  is  used  compared  to  the 
conventional  proton  exchange  membrane. 

2.  Formulation 

The  general  structure  of  an  AEM-DEFC  is  very  similar  to  a  DMFC. 
It  composes  of  a  diffusion  layer  (DL)  and  a  catalyst  layer  (CL)  in  both 
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Cathode  Oxidant 


Fujiwara  et  al.  [6],  2008  Ta 
Hou  et  al.  [7], 2008  Pb 

Modestov  et  al.  [8], 2009  P 

Li  et  al.  [9],  2009  T 

Li  et  al.  [10],  2009  T 

Bianchini  et  al.  [  1 1  ],  2009  T 
Bambagioni  at  al.  [12],  2009  T 
Shen  et  al.  [13],  2010  T 


Pt/Ru 

Pt/Ru/C 

Ru-V/C 


TMPhP-Co/C  100 


Pd— (Ni— Zn)/C  Fe-CO  K 


Zhianietal.  [15],  2011 

Mazin  et  al.  [16],  2011 
Shen  etal.  [17],  2011  1 

Passive 

Bianchini  et  al.  [1 1  ],  2009  T 
Bambagioni  at  al.  [12],  2009  T 
Tadanaga  at  al.  [18],  2010  H1 


FF 


Pd— (Ni— Zn)/C  Fe-CO  K 


J  Tokuyama. 
b  PBI/KOH. 
c  K14  Hypermec™. 
d  Weight  percentage. 
e  Fumasep  FAA. 
f  Natural  (Hydrotalcite, 


U-CO32)  clay. 


anode  and  cathode.  A  hydroxyl  ion  conductive  membrane,  AEM, 
separates  the  anode  and  cathode  compartments.  As  depicted  in 
Fig.  1,  ethanol  diffuses  through  the  anode  DL  to  reach  the  CL.  Using 
the  state-of-the-art  catalyst,  the  full  cleavage  of  a  C-C  bond  in 
ethanol  is  nearly  impossible.  Many  researchers  reported  that  with 
the  available  catalyst,  ethanol  is  only  oxidized  to  acetic  acid  rather 
than  CO2  [6,21],  This  is  advantageous  in  terms  of  reducing  the 
available  CO2  in  the  system  and  mitigating  the  carbonate  precipi¬ 
tation,  though  an  incomplete  EOR  lowers  the  system  energy 
density.  The  EOR  is: 


CH3CH2OH  +  40H  — >CH3COOH  +  3H20  +  4e 


(1) 


While  water  is  produced  in  the  anode,  it  is  required  for  the 
oxygen  reduction  reaction  (ORR)  in  the  cathode: 


02  +  2H20  +  4e~  — >40H-  (2) 

Electrons  of  Eq.  (1)  go  through  an  external  load  and  reach  the 
cathode  DL  and  CL.  The  hydroxyl  ions  of  Eq.  (2)  pass  through  the 
AEM  to  reach  the  anode  CL.  The  overall  reaction  is: 

CH3CH20H  +  02^CH3C00H  +  H20  (3) 

There  are  two  mechanisms  for  ethanol  transport  through  the 
membrane  (crossover):  (i)  ethanol  transfers  through  the  membrane 
from  the  anode  to  cathode  by  the  concentration  gradient,  the  so- 
called  diffusion,  and  (if)  ethanol  also  sticks  to  the  hydroxyl  ions 
(OH- )  traveling  through  the  membrane  from  the  cathode  CL  to  the 
anode  CL,  the  so-called  electroosmosis.  Electroosmotic  coefficient 
refers  to  the  number  of  moles  of  ethanol  attached  to  1  mol  of 
traveling  OH-  through  the  membrane. 

In  the  remainder  of  this  section,  first,  an  analytical  solution  for 
the  electrochemistry  of  an  AEM-DEFC  is  provided.  The  closed  form 
solution  for  two  limiting  cases  of  current  density  is  obtained.  The 
multi-layer  membrane  model  is  then  presented  to  solve  the  mass 
transport.  The  following  assumptions  are  made: 

•  Ethanol  migrated  from  the  anode  to  the  cathode  through  the 
membrane  will  be  completely  oxidized  in  the  cathode  CL.  Thus, 
the  concentration  of  ethanol  in  the  cathode  is  assumed  zero. 

•  The  effect  of  KOH  on  the  kinetics  of  EOR  as  well  as  on  the 
membrane  ion  conductivity  is  neglected. 

•  The  membrane  is  impermeable  to  the  gas  phase. 

•  The  electron  conductivity  of  the  carbon  phase  is  infinity  rela¬ 
tive  to  the  ion  conductivity  of  the  membrane  phase. 

•  Diffusion  and  electroosmosis  are  two  mechanisms  of  ethanol 
transport  through  the  membrane. 

•  Water  concentration  in  the  anode  is  much  higher  than  that  of 
ethanol  and  KOH. 

•  The  kinetics  of  ORR  and  EOR  are  approximated  by  Tafel 
expressions. 
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•  The  concentration  of  ethanol  and  KOH  in  the  anode  channel 
and  the  concentration  of  oxidizer  (O2)  in  the  cathode  channel 
are  constant. 

•  Two  phase  flow  in  the  cathode  is  neglected. 

•  Convective  mass  transport  in  the  DL  and  CL  is  negligible. 

•  The  flow  is  isothermal  and  one-dimensional.  The  change  in 
transport  variables  is  only  considered  in  the  x-direction. 


2,1.  Electrochemistry 

In  the  following  section,  an  analytical  solution  for  the  anode  and 
cathode  overpotentials  is  presented.  The  cell  voltage  for  a  given 
current  density  is  then  obtained.  The  procedure  is  similar  to  what 
was  already  provided  for  a  passive  DMFC  [22,23],  The  EOR  rate  (Rox, 
A  m-3)  in  the  anode  CL  is  estimated  by  Tafel  kinetics: 

)”Kp(w)  (4) 

where  subscripts  a,  e,  and  ox  denote  anode,  ethanol  and  oxidation, 
respectively,  and  superscripts  CL  and  ref  pertain  to  the  catalyst 
layer  and  reference  value.  Jraej  is  the  anode  exchange  current  density 
(A  m3 ),  a  is  the  transfer  coefficient,  77  is  overpotential  (V),  F  is  the 
Faraday’s  constant  (96,485.3,  A  s  mol-1 ),  Ru  is  the  universal  gas 
constant  (8.314,  J  mol  -1K-1),  C  is  the  concentration  (mol  m-3),  T 
is  the  temperature  (K),  and  ya  determines  whether  the  oxidation 

Anode 

a  1 - * - 


reaction  rate  is  zero  or  first  order  with  respect  to  the  ethanol 
concentration  in  the  anode  CL  (C£^): 


Tq  = 


0  cg>crei 
1  ceCLa  <  c?l 


The  anode  overpotential  is: 


(5) 


Va  =  <t>car,a  ~  <t>mem  ~  Ue  (6) 

where  0  is  the  phase  potential  (V),  U  is  the  equilibrium  potential 
(V),  and  subscripts  car  and  mem  denote  carbon  and  membrane 
phases,  respectively.  Note  that  either  the  electronic,  /e(x'),  or 
anionic,  /oh(x'),  current  density  may  be  employed  to  calculate 
overpotentials  at  the  anode  CL.  However,  since  infinite  electron 
conductivity  for  the  carbon  phase  is  assumed  here,  all  the  calcu¬ 
lations  are  based  on  the  anionic  current  density. 

The  anode  CL  depicted  in  Fig.  2(a)  is  the  solution  domain  for  the 
anode  overpotential.  Note  that  the  x'  -coordinate  shown  in  this 
figure  is  different  than  x  -coordinate  in  Fig.  1  which  is  for  the  mass 
transport  calculation.  The  relations  for  the  current  density 
( I0h ,  A  m  2)  inside  the  anode  CL  are: 


0<x!  <  La  a 

I  ,  d<t>mem 

K  ‘OH  =  -ft mem  jgp 


Fig.  2.  A  representative  schematic  of  (a)  anionic,  70 h(x'),  and  electronic,  /e(X),  current  densities  in  the  anode  and  membrane,  and  (b)  single-layer  and  multi-layer  membrane  models. 
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where  a  is  anion  conductivity  of  the  membrane  phase  inside  the 
anode  CL  (CWm-1)- 

The  pertinent  dimensionless  parameters  are: 


where  subscript  CL  pertains  to  catalyst  layer  and  the  stared  values 


<p(W,  o) 


1  +  W,  o 


Attention  is  now  directed  to  the  cathode  overpotential.  Similar 
to  the  anode  side,  the  ORR  rate  in  the  cathode  is  approximated  by 
Tafel  kinetics: 


R  _  ,ref  (Co\  .  c 

Rred  -  J0fl2  I  re/ 


where  subscripts  red,  02  and  c  denote  reduction,  oxygen  and 
cathode,  respectively.  The  variable  yc  determines  whether  the  ORR 
rate  is  zero  or  first  order  with  respect  to  c: 


Upon  combining  Eqs.  (4)  and  (6)  into  Eq.  (7)  and  utilizing  the 
dimensionless  parameters  of  Eqs.  (8)  and  (9),  the  current  density 
inside  the  anode  CL  becomes: 

f!m  +  d^L=  o  (1Q) 

dx>2  dx> 

Eq.  (10)  is  solved  in  the  following  conservative  form: 

dx^  +  4^*0  =  -Rox(x')  +  4(*>  =  — W,o  +  4,o  =  cons 

(11) 

where  subscript  0  denotes  x!  =  0  in  Fig.  2(a)  and  IOH  O  is  the  cell 
current  density.  The  constant  value  in  Eq.  (11 )  must  be  negative  for 
a  physical  solution.  Eq.  (11)  is  solved  for  WW-  Upon  using  the 
condition  of  W(x')  —  0  at  x'  =  1,  the  dimensionless  cell  current 
density  (at  x!  =  0)  is  expressed  as: 

lotto  *-^Rox,o  -  foHfl-tan^Roxv-  W,o^  (12) 

There  is  no  closed-form  solution  for  I0h,o  in  Eq.  (12).  However, 
the  solution  is  obtained  for  two  limiting  cases: 

•1oh,o<s1;  the  tight  hand  side  of  Eq.  (12)  is  small.  For 
0  <  1 8-tan(/3)  <  1 ,  it  may  be  concluded  that  tan(fi)  «  l.Eq.(12) 
is  then  approximated  as  Ioh,o~R<>x,o  -  4,o-  Neglecting  IqH  0 
compared  to  W.o.  the  final  expression  is  further  simplified  as 
1oh,o  =  /W  Considering  the  typical  conditions  of  T =300  K, 
Vmem~3  LaCi~  0.02  mm,  and  a0  =  l,  the  I*0Ha  will  be 

=  7755  A  m-2.  Thus,  the  above-noted  approximation  is  valid 
for  a  wide  range  of  practical  current  density  in  DEFCs,  that  is 
0^1oh,o^700  A  m  2. 

•  W,o  1 :  it  is  assumed  that  due  to  a  very  large  value  of  W, o. 
(Rox, o  -  4,oW71'2/4-  Neglecting  tt2/4  compared  to  fOHO,  Eq. 
(12)  is  simplified  to  I0h, o  =  (Rox.o)°'5-  Followed  by  the  preceding 
discussion,  utilizing  Eqs.  (4)  and  (8),  the  anode  overpotential  is 
obtained  as: 

f  LnQoHfi)  ~  Ln(R0X/R*a)  Wo  «  1 

VaO~<  -  V-  A  -  (13) 

\  2Ln(W,o)  -Ln(R0x/Ra)  W,o»  1 

where  Rox  is  R0x  of  Eq.  (4)  when  rja  —  0.  Eq.  (13)  is  roughly  gener¬ 
alized  for  the  whole  range  of  W,o,  as  follows  [23]: 

Vafl-vQoHfl)  -Ln(W,o)  -  Cn(kox/R*a)  (14) 


coL>coL 


The  cathode  overpotential  is: 

Vc  =  </>car,c  -  -  Uo2  (18) 

Taking  a  similar  procedure  as  Eqs.  (7)-(13),  the  cathode  over¬ 
potential  is  calculated  as: 

(19) 

where  Rred  is  RredofEq.  (16)  when  ?jc  =  0.  Note  that  for  the  ORR,  the 
origin  of  the  coordinate,  x"  =  0,  starts  at  the  interface  of  the 
membrane  and  cathode  CL  towards  the  interface  of  the  cathode  CL 
and  DL.  The  pertinent  dimensionless  parameters  for  the  cathode  are: 


7Ic,0~(pQoH,o)  LnQoHfl )  -  Ln(Rred/R*c ) 


_  W,o 


Rred 


'OH,c 


acFLcCL 


RUT 

=  OcF’ 


K  = 


Lc.CL 


(21) 


In  the  remaining  sections,  for  the  sake  of  simplicity,  Ioh, o  is 
replaced  by  I  denoting  the  cell  current  density. 

The  effect  of  ethanol  crossover  (h"r  e ,  mol  •  m-2  s-1 )  on  the  cathode 
overpotential  is  also  considered.  The  preliminary  crossover  predic¬ 
tions  show  that  parasitic  current  density,  Ip  =  4F-h'^e,  is  small 
compared  to  the  cell  current  density  in  AEM-DEFCs  with  Ce,res  <  3M. 
Assuming  no  mass  transport  limit  for  oxygen  and  ethanol  in  the 
cathode  CL,  Butler-Volmer  relation  can  be  linearized  as  follows  [24]: 


_  Ip/Fc,a  _  Ip/Lc,a 

P,e’C“j0-4F’  P’°2’C_i^-4F 

vf  =  Vc  +  Vi 


(22a) 


+-  ’Ip.O 2,C 


(22b) 

where  yp,e,c  and  VP,Oi,c  denote  the  cathode  overpotentials  due  to  the 
reacted  ethanol  in  the  cathode  CL  and  the  corresponding  oxygen 
reduction,  respectively.  Cell  voltage  (Vceii),  by  definition,  is  the 
difference  between  the  carbon  phase  potential  of  the  anode  (<pCar,a) 
and  the  cathode  (0rarc): 


Vcell  =  0car,c  ~  <t>ca i 


(23) 


Anions  flowing  through  the  membrane  encounter  internal 
resistance  ( Rr,mem )  which  causes  a  potential  drop  in  membrane 
phase  potential  (0mem)  : 

R  ^ mem,q  ^mem.c  _  ^AEM  (24) 

I  Gmem 

Cell  potential  is  presented  as  [22]: 
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B  =  (U02  ~  Ue)  ~Va~  V?  ~  1 1 


(Rcon+LA™) 

\  OmemJ 


'  4F  Lc  CL  Lc  CL 


(33) 


where  RCOn  is  the  contact  resistance  and  Erev  -  (U0l  -  Ue)  repre-  The  effective  diffusivity  coefficient  is  calculated  accounting  for 

sents  the  reversible  cell  voltage.  Erev  is  calculated  as:  the  porosity,  e,  and  the  liquid  saturation,  s,  as: 


AG  AH-TAS 
trev  ~  4f  ~  4 f 


(26) 


where  A  G  is  the  change  in  the  Gibbs  free  energy  of  the  overall 
reaction  of  Eq.  (3).  At  standard  conditions,  Erev  -  1.167  V.  Van’t 
Hoff  relation  is  employed  to  find  the  equilibrium  potential  as 
a  function  of  temperature.  The  change  in  Gibbs  energy  when  the 
temperature  changes  from  T0  to  T  is  [25]: 


AG 

T 


AG0  .  „  (l  1\ 

X+  A 


(27) 


Assuming  T0  =  298  K,  the  equilibrium  potential  as  a  function  of 
temperature  is: 


Erev  =  1.167-3.59  x  lCT4(r-  298) 


(28) 


2.2.  Mass  transport 

Eqs.  (14)  and  (19)  require  ethanol  and  oxygen  concentrations  in 
the  CLs.  Mass  transport  in  a  DEFC  is  similar  to  that  in  a  DMFC  except 
that  the  electroosmotic  transport  in  an  AEM  is  from  the  cathode  to 
anode.  Several  one-dimensional,  analytical  models  for  the  mass 
transport  in  a  DMFC  are  reported,  all  of  which  considered  the 
membrane  as  a  single  layer.  Representative  examples  are 
[22,26,27],  As  will  be  evident,  a  single-layer  model  for  the  AEM 
results  is  a  non-physical  prediction  of  ethanol  transport  through 
the  membrane. 

In  order  to  provide  an  appropriate  model  for  the  ethanol 
transport  through  the  membrane,  a  multi-layer  membrane  model 
is  proposed  in  this  study.  Fig.  2(b)  illustrates  the  difference  between 
a  single-layer  and  a  multi-layer  membrane  model.  Variables  fi'^  e 
and  n'Ee  denote  the  diffusive  and  electroosmotic  flux  of  ethanol, 
respectively: 


%  =  nd,ef  (30) 

where  dCe/dx  is  the  gradient  of  ethanol  concentration  within  each 
representative  layer  (resistance)  shown  in  Fig.  2(b)  and  n de  is  the 
electroosmotic  drag  coefficient  of  ethanol: 


tld,e  =  ndxe  (31) 

where  nd  is  the  electroosmotic  drag  coefficient  of  liquid  within  the 
membrane  and  xe  denotes  the  local  ethanol  mole  fraction.  The 
electrochemical  mass  consumption  of  ethanol  in  the  anode  CL, 
n'ei-che,  is: 


neff  _  f  De-(e-Sa)1-5  anode  CL,  DL 
e  l  De  mem '  emem  AEM  ^ 

Df2  =  D02-(e-(l— Sc))1'5  (35) 


3.  Solution  methodology  and  properties 

Solving  Eqs.  (29)— (33)  provides  the  local  distribution  of  ethanol 
and  oxygen  concentrations.  For  n  layers  (resistances)  in  the 
membrane,  m  in  the  anode  CL  and  k  in  the  anode  DL,  n  +  m  +  k  -  1 
mass  conservation  equations  form  a  well-posed  description  of 
ethanol  concentration  in  the  anode.  The  electrochemical  ethanol 
consumption  in  the  anode  CL  is  uniformly  distributed  to  the  all  m 
layers  in  the  anode  CL  An  iterative  procedure  is  required  to  solve  the 
linear  algebraic  equation  since  the  electroosmotic  drag  of  ethanol 
depends  on  the  ethanol  molar  fraction  (Eq.  (31)).  ATDMA  solver  is 
developed  to  solve  the  linear  equations.  The  iteration  for  the  final 
converged  results  is  checked  for  the  maximum  local  error  criterion: 

max|i/j+1  -  <^10  12  (36) 

where  i  denotes  the  iteration  number,  and  \j/j  is  the  flow  variable  at 
cell  j. 

Distribution  of  oxygen  concentration  in  the  cathode  is  obtained  in 
a  similar  fashion  as  the  ethanol  in  the  anode.  Pure  oxygen  and 
a  constant  concentration  of  ethanol  in  the  cathode  and  anode  channels 
are  respectively  assumed.  The  ethanol  concentration  in  the  cathode  is 
zero.  Membrane  is  supposed  to  be  impermeable  to  the  oxygen.  The 
converged  solutions  of  oxygen  and  methanol  concentrations  are  used 
to  obtain  the  cell  voltage  using  Eqs.  (14),  (19)  and  (25). 

As  will  be  evident,  for  nAEM>  16,  the  ethanol  crossover  predic¬ 
tion  becomes  independent  of  the  number  of  layers  in  the 
membrane.  Unless  otherwise  noted,  all  the  predictions  in  this  study 
are  for  nAEM  =  32,  na  Cl  =  5,  na  DL  =  2,  nc  a  =  5  and  ncDL  =  2. 

Properties  of  AMEs  are  not  well  established  yet,  due  in  part  to  the 
diversity  of  utilized  AEMs  for  DEFC  applications.  There  are  few  studies 
available  regarding  the  measurement  of  transport  properties  in  AEMs 
including  the  diffusivity  coefficient  of  ethanol  through  the  membrane 
[7,19,28,29]  and  the  electroosmotic  drag  coefficient  [19,30].  Table  2 
presents  the  measured  diffusivity  coefficient  of  ethanol  ranging 
from  8.6  x  10-12  to  6.5  x  10-11  m2s_1  for  various  AEMs.  Note  that 
the  table  emphasizes  another  distinct  advantage  of  AEM-DEFCs  over 
DMFCs  in  two  respects:  (i)  ethanol  diffuses  less  than  methanol  in  any 
transport  medium  due  to  its  larger  molecular  structure,  and  (ii)  nearly 
all  the  common  AEMs  exhibit  lower  diffusivity  coefficient  than  Nafion 
no  matter  what  kind  of  fuel  is  utilized  [28],  The  physicochemical 
properties,  cell  geometries  and  operating  conditions  utilized  in  this 
study  are  provided  in  Tables  3  and  4,  respectively. 


n'el-cKe 


4  F-Laja 


(32)  Table  2 

Diffusivity  coefficient  of  ethanol  in  various  AEMs. 


Oxygen  mass  transport  in  the  cathode  is  similarly  solved,  noting 
that  the  membrane  is  impermeable  to  the  oxygen.  The  electro¬ 
chemical  oxygen  consumption  in  the  cathode  CL  n"l_ch02,  accounts 
for  both  the  electrons  from  the  external  load  and  the  crossover  flux 
of  ethanol  through  the  membrane,  n"r  e : 


Type  of  AEM  De  mem(m2s_1)  Investigator 

Tokuyama  2  x  10  11  Varcoe  et  al.  [28] 

PBI/KOH  6.5  X  10  11  Hou  et  al.  [7] 

8.6  x  10~12  Leykin  et  al.  [19] 

PVA/Ti02  composite  polymer  2.81  x  10-11  Yang  et  al.  [29] 
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(-5.7939x1  + 
+1.3679)  x  1 


Diffusivity  in  the  membrane,  ethanol 
Diffusivities,  oxygen 
Ion  conductivity  of  the  membrane 
Electro-osmotic  drag  coefficient 
Anode  transfer  coefficient 
Cathode  transfer  coefficient 
Anode  exchange  current  density 
Cathode  exchange  current  density 
Anode  reference  co 


Cathode  reference  concentration 
Liquid  water  density 
Liquid  ethanol  density 
Contact  Resistance 
Absolute  entropy  (1  atm,  298  K) 


Enthalpy  of  formation  (1  atm,  298  K) 


4 

§ 


1.775  x  10“5(T/273.15)1s: 
3  +  2.57*(T  —  303)  10~2 
2.3  +  0.13 -(T  -303) 


3  x  11T5 

158.0 

69.95 

160.7 

205.33 

-483500 

-285830 

-277690 


Measured  for  a  DMFC 


4.  Results 

The  proposed  model  was  calibrated,  to  the  extent  possible,  by 
predicting  the  experimentally  measured  values  of  the  cell  voltage 
and  power  density  [10],  Fig.  3  presents  a  comparison  of  the  cell 
voltage  solution  of  Eq.  (25)  to  the  polarization  curve  and  the  power 
density  of  [10]  when  Ce,res  =  1  M  and  pure  oxygen  is  fed  to  the 
cathode  channel.  The  cell  temperature  is  T  =  303  K.  A  good 
agreement  is  observed  between  the  numerical  predictions  and  the 
experimental  data. 


4.1.  Electroosmotic  drag  coefficients 

The  focus  of  this  section  is  to  emphasize  the  differences  in  mass 
transport  within  the  membrane  corresponding  to  the  positive  and 
negative  electroosmotic  drag  coefficients.  Electroosmotic  transport 
of  one  species  in  the  membrane  linearly  scales  with  the  local  mole 
fraction  of  that  species,  Eqs.  (30)  and  (31).  In  most  DMFC  models 
developed  so  far,  as  a  typical  system  in  which  the  electroosmotic 
drag  coefficient  is  positive,  the  membrane  is  considered  a  single 
layer  and  the  electroosmotic  transport  of  methanol  is  modeled 
based  on  the  concentration  of  methanol  at  the  interface  of  the 
anode  CL  and  membrane  [22,36—42],  This  approach  works  properly 
for  DMFCs  where  the  directions  of  both  electroosmotic  and 


id  operating  parameters. 


Parameters 


Anode  DL  thickness 
Anode  CL  thickness 

Cathode  CL  thickness 
Cathode  DL  thickness 
Porosity  of  DL  and  CL 
Porosity  of  membrane 
Cell  temperature 
Anode  liquid  saturation 
Cathode  liquid  saturatior 


I  (A  nr2) 


membrane  model  ( n MM  =  2,  4,  8,  16). 
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Fig.  5.  The  effect  of  positive  and  negative  electroosmotic  drag  coefficients  on  the 
various  ethanol  crossover  mechanisms  as  well  as  the  total  crossover,  0  <1  < 
1000  Am“2. 


diffusive  transport  are  from  the  anode  to  the  cathode.  However,  for 
a  fuel  cell  system  employing  an  AEM  where  the  electroosmotic 
transport  is  negative  (from  the  cathode  to  anode)  following  the 
similar  procedure  fails  and  results  in  a  non-physical  prediction  of 
ethanol  crossover. 

In  an  AEM-DEFC,  the  electroosmotic  mechanism  operates  in  the 
opposite  direction  of  the  diffusion.  The  net  mass  flux  of  ethanol  is 


o  too  200  300 


ethanol  distribution  in  the  anode  and  membrane,  and  (b)  on  the  various  ethanol 
crossover  mechanisms  within  the  membrane,  when  Ce,res  =  1  M  and  /  =  500  A  m-2. 


Fig.  7.  Diffusive  (dashed  line),  electroosmotic  (dashed-dot  line),  and  total  (continuous 
line)  ethanol  crossover  through  an  AEM,  when  Ce-res  is  0.5, 1,  2,  and  3  M. 


therefore  obtained  by  subtracting  the  electroosmotic  transport 
from  the  diffusive  transport  (Fig.  2(b)).  An  appropriate  number  of 
layers  in  the  membrane  is  required  to  get  a  physical  predictions  of 
the  ethanol  transport  through  the  membrane.  Since  the  ethanol  is 


o  too  200  300 

x(|im) 


b  Cathode 


Fig.  8.  Local  distribution  of  (a)  the  ethanol  concentration  in  the  anode  and  the 
membrane,  and  (b)  the  oxygen  concentration  in  the  cathode,  for  various  values  of 
current  density. 
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Fig.  9.  Local  distribution  of  the  diffusive  (dashed  line),  electroosmotic  (dashed-dot  line),  and  total  (continuous  line)  ethanol  flux  in  the  membrane  of  an  AEM-DEFC,  when  Cei„,s  and  / 
are  respectively:  (a)  1  M  and  100  A  m~2,  (b)  1  M  and  500  A  nr2,  (c)  3  M  and  100  A  nr2,  and  (d)  3  M  and  500  A  nr2. 


supplied  in  the  anode,  the  negative  ethanol  crossover  is  not 
a  physical  solution.  Fig.  4  delineates  the  predictions  of  the  ethanol 
crossover  through  the  membrane  for  various  selected  number  of 
layers  in  the  membrane,  nAEM,  when  naa  =  5,  naDL  =  2  and 
Ce,res  =  1  M.  The  conventional  one-layer  membrane  model,  rep¬ 
resented  by  nAEM  =  1,  results  in  a  negative  ethanol  crossover  at 
medium  and  high  current  densities.  Increasing  the  number  of 
layers  in  the  membrane  mitigates  the  nonphysical,  negative 
ethanol  crossover  and,  as  evident,  the  predictions  become  grid 
independent  for  nAEM  >16. 

Fig.  5  delineates  the  difference  between  various  ethanol  cross¬ 
over  mechanisms,  and  the  total  crossover  in  a  DEFC  for  the  positive 
and  negative  electroosmotic  drag  coefficients  ( nd ).  The  positive  nd 
represents  a  DEFC  using  a  proton  exchange  membrane.  All  condi¬ 
tions  are  identical  for  two  cases  of  positive  and  negative  nd. 
Reported  values  are  averaged.  While  the  difference  in  the  diffusive 
transport  for  two  cases  of  negative  and  positive  nd  is  nearly  negli¬ 
gible,  a  significant  difference  is  observed  in  electroosmotic  trans¬ 
port  and  the  total  crossover.  The  trend  of  predicted  total  crossover 
for  the  positive  nd  case  shows  a  good  agreement  with  the  trend  of 
both  experimental  measurements  [4]  and  numerical  predictions 
[43,44]  of  methanol  crossover  in  a  DMFC,  a  representative  of 
a  system  with  nd  >  0. 

Fig.  6  illustrates  the  influence  of  the  change  in  electroosmotic 
drag  coefficient  on  the  local  distribution  of  ethanol  concentration  in 
the  anode  and  membrane  and  the  local  ethanol  flux  within  the 
membrane,  when  Ce, res  =  1  Mandf  =  500  A  m-2.  Fig.  6(a)  shows 
the  difference  between  the  ethanol  concentration  distribution  in 
the  anode  and  the  membrane  for  positive,  negative  and  zero  elec¬ 
troosmotic  drag  coefficients.  It  is  observed  that  the  greatest 
difference  in  ethanol  concentration  occurs  within  the  membrane. 
While  Fig.  5  demonstrates  that  the  averaged  value  of  diffusive 
transport  of  ethanol  in  the  membrane  is  nearly  identical  for  the 
positive  and  negative  nd,  Fig.  6(a)  delineates  that  the  local  diffusive 
transport  of  ethanol  within  the  membrane  should  significantly  be 


different  for  various  values  of  nd.  For  nd  >  0,  the  gradient  of  ethanol 
concentration,  dCe/dx,  varies  from  a  near  zero  value  at  the  interface 
of  CL  and  AEM  to  a  large  value  at  the  interface  of  the  AEM  and 
cathode.  When  nd  <  0,  dCe/dx  within  the  membrane  decreases  as  x 
increases.  Fig.  6(b)  depicts  the  local  diffusive,  electroosmotic  and 
total  ethanol  flux  within  the  membrane  for  the  positive  and 
negative  nd.  As  evident,  the  diffusive  flux  of  ethanol  for  nd  >  0  is 
maximum  where  it  is  minimum  for  nd  <  0,  and  vice  versa.  The 
magnitude  of  electroosmotic  fluxes  for  both  positive  and  negative 
nd  nearly  follow  the  same  trend,  and  change  from  a  maximum  value 
at  x  =  280  pm  to  zero  at  x  —  310  pm.  The  electroosmotic  drag 
within  the  membrane  is  a  function  of  ethanol  concentration.  A 
significant  difference  between  the  net  flux  of  ethanol  in  the 
membrane  is  observed  for  two  cases  of  positive  and  negative  nd. 
The  net  ethanol  flux  is  132.7  x  10-6  mol  m-2s-1  for  nd  —  2.3, 
whereas  it  drops  to  7.5  x  10-6  mol  m_2s_1  for  nd  =  -2.3. 


4.2.  Mass  transport  in  AEM-DEFCs 

Attention  is  now  directed  to  the  AEM-DEFCs  where  nd  <  0.  Fig.  7 
displays  averaged  values  for  diffusive,  electroosmotic  and  total 
ethanol  crossover  through  the  membrane,  when  0.5  <  Ce  res  <  3  M. 
Diffusive  crossover  of  ethanol  monotonically  decreases  with 
increasing  current  density.  It  is  increasing  for  higher  concentration 
of  ethanol  in  the  anode  reservoir.  Electroosmotic  crossover  of 
ethanol  from  the  cathode  to  anode  is  zero  where  /  =  0.  It  increases 
to  a  maximum  value  due  to  the  increasing  current  density  and  then 
starts  to  decrease  to  a  minimum  at  f=4m  owing  to  the  low 
concentration  of  ethanol.  Electroosmotic  crossover  is  also 
increasing  as  more  concentrated  ethanol  is  fed  to  the  anode.  For  all 
four  cases  shown  in  Fig.  7,  the  net  ethanol  crossover  significantly 
decreases  from  the  maximum  value  at  /  =  0  to  a  near  zero  value  at 
1  >  500  A  m-2 .  Reduction  in  ethanol  crossover  is  the  unique  feature 
of  an  AEM-DEFC. 
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Fig.  8  delineates  the  local  distribution  of  ethanol  in  the  anode 
and  membrane  as  well  as  the  distribution  of  oxygen  in  the  cathode, 
for  100  <  /  <  900  A  m-2  and  Ceres  —  1M.  As  expected,  both  ethanol 
and  oxygen  concentrations  decrease  as  current  density  increases. 
An  initial  increase  and  then  a  decrease  in  the  electroosmotic  drag  of 
ethanol  from  the  cathode  to  anode  in  Fig.  7  is  adequately  explained 
through  the  ethanol  distribution  in  the  membrane.  Electroosmotic 
drag  is  zero  when  1  =  0  and  increases  for  increasing  I  (Eq.  (30)). 
However,  Ce  in  the  membrane  reduces  as  I  increases,  which  in  turn 
decreases  the  electroosmotic  drag  coefficient,  according  to  Eq.  (31 ). 
Electroosmotic  drag  from  the  cathode  to  the  anode  is  minimum  at 
I -him  due  to  the  depletion  of  membrane  from  ethanol  (Fig.  8(a)). 
Compared  to  Fig.  8(a),  Fig.  8(b)  depicts  that  the  oxygen  concen¬ 
tration  in  the  cathode  never  falls  below  the  range  that  causes  mass 
transport  limitations. 

Fig.  9  shows  the  local  electroosmotic  and  diffusive  fluxes  of 
ethanol  within  the  membrane,  when  I  =  100,  500 Am-2  and 
Ce,res  =  1 ,  3  M.  Electroosmotic  flux  of  ethanol  from  the  cathode  to 
anode  is  maximum  at  x  =  280  pm,  where  the  Ce  has  the  highest 
value  in  the  membrane.  It  is  zero  at  x  =  310  pm,  where  it  is 
assumed  that  the  concentration  of  ethanol  is  zero.  The  local  diffu¬ 
sive  transport  of  ethanol  also  starts  from  a  maximum  at 
x  =  280  pm,  and  ends  at  a  value  equal  to  the  net  total  flux  of 
ethanol  at  x  =  310  pm.  As  evident,  the  local  magnitude  of  elec¬ 
troosmotic  flux  of  ethanol  in  the  membrane  roughly  scales  with  the 
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Fig.  11.  The  effect  of  increase  in  the  membrane  thickness  on  the  local  distribution  of  (a) 
the  ethanol  concentration,  (b)  the  electroosmotic,  and  (c)  the  diffusive  flux  of  ethanol 
within  the  membrane,  when  CE  rH  =  1  M,  1  =  500  A  m-2  and  0.03  <  X^em  <  0  09  nun. 

magnitude  of  the  diffusive  transport  of  that,  and  the  trend  for  both 
transport  mechanisms  is  nearly  the  same. 

Poor  ion  conductivity  of  the  available  AEMs  is  one  of  the  major 
concerns  in  the  development  of  AEM-DEFCs.  Considering  the 
significant  decrease  in  the  ethanol  crossover,  especially  at  the 
higher  current  densities  (Fig.  7),  decreasing  the  thickness  of  the 
membrane  is  one  of  the  feasible  and  effective  solutions  to  the 
above-noted  problem.  Fig.  10(a)  shows  the  influence  of  membrane 
thickness  on  the  ethanol  crossover  when  Ce,res  =  1  M  and  0.03  < 
Laem  <  0.09  mm.  It  is  illustrated  that  reducing  the  thickness  of  the 
membrane  results  in  a  significant  increase  in  ethanol  crossover 
when  f-»0.  Fig.  10(b)  shows  the  effect  of  membrane  thickness  on 
the  cell  polarization.  At  low  current  densities  (/-» 0),  the  cell 
performance  is  deteriorated  with  decreasing  membrane  thickness. 
It  is  attributed  to  the  higher  cathode  overpotential  due  to  the 
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increased  ethanol  crossover  for  the  thinner  membranes.  The 
performance  of  the  cell  with  thinner  membrane,  however, 
increases  at  the  higher  current  densities  for  two  reasons:  (i)  lower 
resistance  against  OH~  transport  through  the  membrane  (Eq.  (24)), 
and  (ii)  negligible  ethanol  crossover.  Therefore,  reducing  the 
thickness  of  the  membrane  is  an  effective  method  to  compensate 
for  the  poor  ion  conductivity  of  AEMs,  if  the  DEFC  operates  at  high 
current  densities. 

Fig.  11  shows  the  local  distribution  of  ethanol  concentration, 
electroosmotic  drag  and  diffusive  transport  of  ethanol  in  the 
membrane  for  Ce,res  =  1  M,  /  =  500  A  rrr2  and  0.03  <  LAEM  < 
0.09  mm.  The  total  ethanol  flux,  h",  made  up  of  the  convective 
(electroosmotic)  and  diffusive  fluxes  is: 

n"  =  Ceii  -  Df  ^  (37a) 

“  ~%(Clko~  Ce  .  CK0„)  (37b) 

Peclet  number,  Pe,  is  defined  as  the  ratio  of  the  convective  to  the 
diffusive  flux: 

Pe  =  .  uLaem  (38) 

^e,mem  ’  £mem 

As  evident  in  Fig.  11(a),  ethanol  flow  through  the  membrane 
becomes  diffusion  dominant,  low  Pe,  as  the  membrane  thickness 
decreases.  Fig.  11(b)  and  (c)  illustrate  that  both  electroosmotic  and 
diffusive  drags  of  ethanol  in  the  membrane  increase  upon 
a  decrease  in  the  membrane  thickness.  The  electroosmotic  trans¬ 
port  and  Ce  distribution  in  the  membrane  follow  the  similar  trend, 
while  diffusive  transport  scales  with  the  local  derivative  of  Ce. 


5.  Conclusions 

A  one-dimensional,  isothermal,  multi-layer  membrane  model  is 
proposed  to  study  the  mass  transport  in  a  DEFC  using  an  AEM.  The 
electrochemistry  is  solved  analytically  by  assuming  Tafel  expres¬ 
sions  for  both  EOR  and  ORR.  Explicit  expressions  for  the  anode  and 
cathode  overpotentials  are  obtained  for  two  limiting  cases.  From 
the  proposed  model,  the  following  conclusions  are  made: 

•  Employing  a  multi-layer  membrane  model  is  essential  to 
provide  a  physical  and  proper  description  of  mass  transport  in 
an  AEM-DEFC. 

•  A  substantial  difference  in  ethanol  distribution  within  the 
membrane  exists  for  the  positive  and  negative  electroosmotic 
drag  coefficients. 

•  A  significant  decrease  in  ethanol  crossover  is  observed  upon 
using  an  AEM  compared  to  using  a  proton  exchange 
membrane. 


Reducing  the  membrane  thickness  is  an  effective  and  reliable 
method  to  compensate  for  the  low  ion  conductivity  of  AEMs,  if 
the  DEFC  operates  at  high  current  densities,  f>500  A  m~2. 
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